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Abstract. The spline space C^(A) attached to a subdivided do- 
main A of W'' is the vector space of functions of class C" which 
are polynomials of degree < fc on each piece of this subdivision. 
Classical splines on planar rectangular grids play an important 
role in Computer Aided Geometric Design, and spline spaces over 
arbitrary subdivisions of planar domains are now considered for 
isogeometric analysis applications. We address the problem of de- 
termining the dimension of the space of bivariate splines C^(A) 
for a triangulated region A in the plane. Using the homological 
introduced by Billcra ( 19881) . we number the vertices and estab- 
lish a formula for an upper bound on the dimension. There is no 
restriction on the ordering and we obtain more accurate approxi- 
mations to the dimension than previous methods and furthermore, 
in certain cases even an exact value can be found. The construction 
makes also possible to get a short proof for the dimension formula 
when k > Ar + I, and the same method we use in this proof yields 
the dimension straightaway for many other cases. 



1. Introduction 



Let A be a connected, finite two dimensional simplicial complex, 
supported on |A| C M^, with |A| homotopy equivalent to a disk. We 
denote by C^(A) the vector space of all functions on A that, re- 
stricted to any simplex in A, are given by polynomials of degree less or 
equal than k. These functions are called splines, and they have many 
practical applications, including the finite e lement method for solving 
differential equations ( Fix and Strang) ( 20081 )). Recently they have also 



been considered for isogeometric analysis applications (jCottrell et al. 
(120091 ) )■ i n corn p uter-a ided design for modeling surfaces of arbitrary 
topology (IFarinl (119931 )). In these application areas, sphne functions 



of low degree are of particular interest, and the degree of smoothness 
attainable is an important design consideration. 
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A fundamental problem is to determine the dimension of this vector 
space as a function of known information about the subdivision A. A 
serious difficulty for solving this problem is that the dimension of the 
space Cl{A) can depend not only on the combinatorics of the subdi- 
vision, but also on t he geometry of the triangulation i.e., how |A| is 
embedded in M^. In Schumaker (jl984), the author presented a lower 
and an upper bound on the dimension of spline spaces of arbitrary de- 
gree and smoothness for general triangulations; using Bernstein-Bezier 
methods, a formula fo r the d imension for k > Ar + 1 was obtained by 



m 



HongI (jl99ll ). and in 



Alfeld and Schumaker! ( 1987D. The result was extend ed to k > 3r + 2 



Alfeld and Schumaker! (!l990l ) a dimension for- 
mula is proved for almost all triangulations for k = 3r + 1. How- 
ever, there are no explicit formulas for the dimension of the spline 
spaces C^(A) for degree < 3r + 2 for general triangulation s. The 
use of homological algebra in spline theory was introduced by !Billera 
( 11988! ). He o btained the dimensioii of Cl(A ) for all k for generic tri- 
angulations. Schenck and StillmanI f !l997a} ). introduced a chain com- 
plex different from that used by Billera; this complex was s t udied by 
Schenck anc Stillman! (!l997a}). !Schenck and Stillma"^ f!l997bf). !Schenck! 



( 11997! ). and !Geramita and Schenck! ( 11998! ). The lower homology mod- 
ules of the chain complex in this construction, differ from the one in- 
troduced by Billera, and they have nicer properties. The connection 
between fat points and the spline space defined on A, allows to give 
a complete characterization of the free resolutions of the ideals gener- 
ated by power of linear forms appearing in the chain complex, and so 
to prove the dimension formul a for C^(A) for sufficiently high degree 
( Schenck and Stillman ( 1997bl )). 

The main contribution of the paper is a new formula for an upper 
bound on the dimension of the spline space. The formula applies to any 
ordering established on the inte rior vertices of the partition, contrarily 
to the upper bound formulas in !Schumaker! (!l984[ ). !Lai and Schumaker 
( 12007! ). Having no restriction on the ordering makes it possible to 
obtain accurate approximation to the dimension and even exact value 
in many cases. As a consequence, we give a simple proof for dimension 
formula when k > Ar + 1. 

The paper is structured as follows. In Section !5] we recall the con- 
struction of this chain complex and some of the properties of the ho- 
mology modules. We describe the dimension as a sum of a combinato- 
rial part and the dimension of an homology module. This latter term 
which happens to be zero in high degree is always non-negative, so the 
formula gives a lower bound for the dimension of the spline space for 
any degree k. We present the main result in Section !3l which is an 
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upper bound on the dimension of the sphne space. In Section H] we 
compare the formulas o f the l ower and the upper bound with those 
appearing in ISchumakeii (jl984l ). The result about the exact dimension 



for /c > 4r + 1 is proved in Section [51 This latter result and some 
other examples that we present in Section [6l illustrate the interest of 
the homology construction for proving exact dimension formulas. 

2. Construction of the chain complex 



We r eproduce some notations and definitions presented in lSchenck and Stillman 
fll997al ). restricting them to the case where A is a planar simplicial 
complex supported on a disk. 

Denote by A° the set of interior faces of A, and by A^ {i = 0,1,2) 
the set of i-dimensional interior faces whose support is not contained 
in the boundary of |A|. We denote by /° the cardinality of these sets, 
and by dA the complex consisting of all 1-faces lying on just one 2- 
face as well as all subsets of them. As it will be convenient to study 
the dimension of the vector space C^(A), we embed A in the plane 
{z = 1} C and form the cone A over A with vertex at the origin. 
Denote by C^(A) the set of splines on A of C''"-smoothness and degree 
exactly k. Then C"'(A) := ©>oC^(A) is a graded M-algebra and there 
is an isomorphism of R-vector spaces between C ^(A) and the elements 
in C"'(A) of degree exactly k f Billera and Rose f 199l[ l). in particular 



dimCUA) = dimCUA). 

Define R := ]R[x, For an edge r G A^, let ir denote a non-zero 
homogeneous linear form vanishing on f, and define the ideal j7'(/3) of 
R for each simplex /3 G A'^ as follows: 

J{cx) = (0) for each a G A° 

J{t) = for each r G A? 

J{-f) = (C')r97 for each 7 g A°, r G A?. 

Consider the chain complex TZ defined on A'' as TZi = R^^ ^ with the 
usual simplicial boundary maps di used to co mpute the rel ative (mod- 
ulo d/S) homology with coefficients in R as in iBilleral f lTgSSh . Let 11/ J 
be the chain complex obtained as the quotient of 71 by J , 

7^ ^ n/j{r) % 7^/ j(7) 











o-eA2 



tga; 



where the maps di are the induced by the simplicial boundary maps di. 
The complex TZ/ J was introduced by ISchenck and StillmanI f!l997ah . 
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and agrees with the com p lex st udied by Billera except at the vertices. 
It was shown in Billeral ( 1988 ) that C"'(A) is isomorphic to the top 
homology module oilZj J 

H2ijl/J) ■= ker(92). 

Let us consider the short exact sequence of complexes 

— y J — yU — y n/J — y 

that gives rise to the long exact sequence of homology modules 

-> H2{n) ^H2{n/j) Hi{j) Hi(n) 

(1) ^ H,{n/j) ^ HoiJ) ^ Hoin) ^ Hoin/j) ^ o 

Since A is supported on a disk, both Hq{1Z) and HiilZ) are zero. Hence, 
by the long exact sequence, it implies that Ho(Jl/J') is also zero and 
Hiiji/J) is isomorphic to Ho{J). Applying the Euler characteristic 
equation to the chain complex 71/ J 

x{H{n/j)) = xWJ), 

and considering the modules in degree exactly k, leads to the formula 

2 

(2) dimCI(A) = ^(-ir 5^ dunn/J{f3)k + dunHo{J)k. 

The aim is to determine the modules in the previous formula as func- 
tions of known information about the subdivision A. 
We know that 

2 



(3) 



(4) 



J2 dimn/J{T)k = 



k + 2 
2 



k + 2-{r + r 
2 



For the computation of dim7?./j7'(7)fc, ISchenck and StillmanI (Il997bl ) 
proposed the following resolution for 7l/J^{'y)k- Without loss of gener- 
ality we translate 7 to the origin and assume that the linear forms in 
J7'(7) involve only variables x,y. Letting ii'^^, . . . be a minimal 
generating set for Jipf), then a free resolution for IZ/ J{pf)k is given by: 

(5) ^ n{-VL - If e n{-Q.f ®]=in{-r - i)^n nijij) o 

where — 1 is the socle degree oilZ/ J{^)] VL and the multiplicities a 
and h are given by 

+1, a = t(r + 1) + (1 b = t-l 



(6) 



t- 1 



a. 
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In the case t = 1, we take A = b = Q = Oso that 7^(-fi - 1)" © 
71{—Qy = 0. Applying this to each vertex 7j G Aj], and letting tj, Qi, 
tti and bi denote the values for t, fl, a and b at 7j respectively, leads to 

fO 

Jo r 



r / 

dim 7^/J(7.)fc = ( 



A; + 2 



U 



k 



2 - (r 
2 



+ 



(7) 



A; + 2- (a + r 
2 



A; + 2 - 
2 

Let us notice that the previous formula does not change when we take 
ti as the number of different slopes of the edges containing the vertex 
7i. Then ©, ©, and ([7D, and the fact that dim Hi{n/ J) k > 0, yields 
the following theorem: 

Theorem 2.1. The dimension o/C^(A) is bounded below by 



dimC[(A)> 



k + 2 



-E 

i=l 



k + 2-{r + r 
2 



/c + 2 - (r + i; 
2 



k + 2 -a, 

2 



k + 2-{ni + r 

2 



where ti is the number of different slopes of the edges containing the 
vertex 7j, and 

tiV 



ni 



ti 



+ 1, 



ti (r + 1) + (1 - ti) a 



t.:-l 



Qi. 



In Schenck fll997[ ). it was proved that the homology module Hi(TZ/ J)k 
vanishes for sufficiently high degree, thus the lower bound in the latter 
theorem is actually the exact d imension formula for C^(A) when A; 3> 
(iGeramita and Schenckl . [l998l . Theorem 4.2). 

3. An upper bound on the dimension of Q(A) 

Let us fix an ordering 71, . . . ,7^0 for the vertices in Aq. For each 
vertex 7,, denote by A^(7j) the set edges that contain this vertex, and 
define ti as the number of different slopes of the edges connecting 7j to 
one of the first i — 1 vertices in the list or to a vertex on the boundary. 

Theorem 3.1. The dimension of Cl{A) is bounded by 



dim CliA) < 



k + 2 



A; + 2 - (r + 1' 



A; + 2 - (r + 1) 
2 

k + 2-ni 
2 



k + 2-{ni + i) 
2 
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with ii as we have defined above and 

+ 1, di = ii (r + 1) + (1 - ii) Qi, bi = ii - 1 - di. 
if ti> 1 and di = bi = Qi = if ti = 1 . 

Proof. By the long exact sequence ([T]), and the fact that HiilZ) = 0, 
we have the short exact sequence 

H^iTZ) H^iJZ/J) H^{J) 

The Euler characteristic equation apphed to the complex TZ implies 
that H2{TZ)k = Rk- Hence the isomorphism between C"'(A) and the 
homology module H2{1Z/ J) implies that 

dimC^(A) = AimRk + <limHi{J)k 

where the complex of ideals J (as defined above) is given by 

^ © J{r) % J{l) 

and where Hi{J) = ker^i = Ki. Define Wi := lm{di). By the exact 
sequence 

i^i ^ J{r) — > 

we get 

dimC^(A) = dimi?fc + ^ dim J(r)fc - dim(iyi)fc 

reA? 

Therefore, to find an upper bound on dimC^(A) it is enough to find 
a lower bound on the dimension of Wi = Imdi in degree k. Let us 
consider the map 

-^i: J{r)[r]^ R[t\^] 

r=(7,7')eA0 7,GA[jreA'(7) 

such that 5i([r]) = [r\'y] — [t\j'] for r = (7,7') G A^, and the map 

^1^0 R[r\i]^®Rb] 

7,gA[jrgAr(7) 7gA0 

with 

y,i([r|7]) = [7] if 7eA°, 
= if7^A°. 



Qi 



Ur 
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Then, we have di = (fi o Si. We consider now the map 
^1^0 ^[^M^ R[r\l] 

with 7ri([r|7]) = if 7 is the end point of biggest index of r and 
7ri([r|7]) = [r|7] otherwise. We denote by di = ipi o tti o Si. 

For 7 G Aq, we denote by A^(7) the set of interior edges r con- 
necting 7 to another vertex which is not of bigger index. Let J'{'y) = 
St67V(7) RK^^ J' in)- By construction, we have 

lm~di= J(7) [7] 
7eA0 

and dim(14^i)fc = dim(Im9i)fc > dim(Im(9i)fc. 

Using formulas ([3]), (jlj) and the resolution ([5]) applied to J'i'j), we 
obtain the upper bound 



dimCI(A)< + 



k + 2\ „n /A; + 2 - (r + 1; 

2 



fO 

E 



with tj = |A^(7i)| and fij, cij, 6, defined as in (Q. 

As an immediate consequence of this theorem we mention the fol- 
lowing result. 

Corollary 3.2. If for a numbering on the vertices in Aq, either ti = ti 
orti > r + 1 for every vertex 'ji, then the upper bound we get, equals the 
lower bound, obtaining so the exact dimension formula for the spline 
space. 

Proof. We compare the terms corresponding to each interior vertex in 
both formulas. If ti = ti then they are trivially the same. By definition 
([6]), for t > r + 1, Q and a are both constant equal to r + 1. Hence 
if r + 1 < ti (< ti), the terms in the binomials are the same when 
we compare them in the formulas in Theorem 12.11 and 13. 1[ Since the 
respective terms in ti and U cancel out then we get the equality of the 
bounds. 
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4. The bounds on dimC^(A) given by ISchumakerI ( 1 19841 ) 

In t his section we compare the bounds on dim C^(A) found by lSchumakei 
( ll984j ). with the lower and upper bounds given in the previous two sec- 
tions. With the notation as before, the upper bound presented by 
Schumaker can be stated as follows. 



Theorem 4.1. nSchumaken . \1984 - Theorem 2.1) Suppose that the ver- 
tices 7, G Aq of the partition are numbered in such a way that each pair 
of consecutive vertices in the list are corners of a common triangle in 
A. For each 7j define ti as the number of edges with different slopes 
joining the vertex 'ji to a vertex in the boundary of A or to one of the 
first i — 1 vertices. Then 



(8) 



Jo 



k + 2 



r + 2 



fo k~ 



In the same article, the author also presented a lower bound on 
dimC^(A). The formula can be obtained by replacing tj by tj in ([8]). 
We shall prove that that formula for the lower bound is the same as the 
formula we presented in Theorem 12.11 Essentially the same proof shows 
that the upper bound formulas, the one in the previous theorem and 
the one we presented in Theorem (13. ip . coincide. With the exception 
that the upper bound presented by Schumaker can only be applied for 
certain numberings on the vertices. That restriction makes sometimes 
not possible to find an upper bound. In Section IH] we will include some 
examples of those situations, and also some cases when not having that 
restriction leads to find the exact dimension of the space. 

Remark. The above cited article defined ti as the number of slopes 
of the edges containing 7^ but not containing any of the i — 1 vertices 
in the list. This coincides with the definition of tj we consider here, 
except that the reverse counting on the vertices in Aq is used. 

Lemma 4.2. Let t >2 and Q defined as in ^ above. The notations 
in the upper bound formulas (and the lower bounds, respectively) are 
linked in the following way: if Q = r + i then r + j + 1— j-t>0 
only when j < i — I. 

Proof. For any value of t > 2 we have 
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Thus, r + l<^7<2r + l, and we can write Q = r + £ for some 
integer £ between 1 and r + 1. The interval of t where the value of Q 
is r + £ is given as follows: 



(9) 



I±£<t<m£_i2 when £>2 



t > r + 1 when £ = 1. 



On the other hand, for fixed r and t, the condition 
(10) r+j + l-j-t>l 

is satisfied if and only if t < (r + ])/]■ Thus, the biggest number j 
subject to condition ( ITO|) must satisfy 

r + {j + 1) ^ ^ ^ r+j 

J+1 ~ 3 ' 

From ([9]) , the previous relation holds if and only if = r + (j ' + 1) . Let 
t > 2 and Q, a and b be defined as in ([6]). We consider the following 
two formulas: 

^^^^ a + 2-{r + i)\ _^rk + 2- _^fk + 2- {n + r 



2 \ 2 \ 2 




Lemma 4.3. If t > r + 1 and k > r then the formulas 01]) anc? / fl^) 
are egwa/. 

Proof. By dH]) we know that if t > r + 1 one has Q = r + 1. From the 
definition in (Q, we have a = r + 1 and 6 = t — (r + 2). Thus, f|TT]) 
reduces to 

{0 if A; + 1 - r < 2 (A; = r) 

r + 2 ifA;-r = l 

('^') - C?) otherwise. 

This can be simplified to C^^^) — ('^^^) for any value of k and r. By 
Lemma 14.21 the relation (JTOl) in this case is not satisfied by any value 
of j, then also ( IT^ simplifies to that expression. 

Lemma 4.4. Let t be an integer >2. If t <r + 1 and k > r then the 
formulas ( fTTj) and ( fi^) are egwa/. 
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Proof. By Lemma 14.21 we know that Q = r + i for some integer 
i > 2. Depending on the value of i the binomials in flTT]) could be 
automatically zero. We need to consider three possible situations k — 
r>i+l, k — r = l and k — r < i. The formula reduces to a different 
expression in each case, but it is not difficult to check that they are 
respectively equivalent to the expressions we get from f|T2|) . 

Proposition 4.5. The lower bound on dimC^(A) give n in Theorem 



\2.1\ coincides with the lower hound in \Schumaken ( fj 984) [ Theorem 3.1 
p. 256]. 

Proof. Since any interior vertex of A is contained in at least two edges 
with different slopes, then ti > 2 for every i = 1, 2, ... /q. Collecting 
the vertices 7^ with the same value ti, the statement follows directly 
applying Lemmas 14.31 and for those values of ti. 

Proposition 4.6. // the vertices in are numbered as in Theorem 
d^. then the upper boun d on dim CUA) giv en in Theorem \3.1\ coin- 



cides with the bound ^ in \Schumaken ^1984) [ Theorem 2.1 p. 252] 



Proof. This follows similarly as in the latter proposition. By collecting 
the vertices in A" with the same value ti, we apply Lemmas 14.31 and 
4.41 with the values ti. The only case that remains to be considered is 
ii = 1. It corresponds to the third term in the formula in Theorem 
flO) . Since 



k + 2 - (r + 1)\ _ fk + 2\ + ^ 



k—r 

2 \ 2 \ 2 ' 

j=i 

the statement follows. 

We summarize the results in the following theorem. Let us denote 
by LBH, UBH and LBS, UBS the respective lower and up per bounds 



obtained by Theorems 12.11 and 13.11 and the ones obtained in lSchumaker 



(119841 ). 



Theorem 4.7. Let A be a connected, finite two dimensional simplicial 
complex, supported on a disk. Then, 

LBH (dim C^- (A)) = LBs(dimC^(A)) 

UBH (dim CI,' (A)) < UBs(dimC^(A)). 

Proof. This follows from Propositions 14.51 and 14. 6[ and the fact that 
the formula in Theorem 13.11 can be applied to any numbering on the 
interior vertices of A. 
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5. Dimension formula for degree A; > 4r + 1 

In this section we present an alternative proof for the dimension 
formula of the spline space C^(A) for k > 4r + 1. The proof is substan - 



tially shorter than the one presented in lAlfeld and Schumakerl (119871 ) . 
We include the following notation: let Aj and Af be (respectively) the 
sets of z-dimensional faces and i-dimensional boundary faces. Denote 
by /i(A) and /f (A) the cardinality of the preceding sets, respectively. 
We begin by stating the following lemma which we will need later on. 



Lemma 5.1. /(Schenck and Stillman . \l997d . Lemma 3.3) If ^ is a 



triangulated region in M^, then there exists a total order on Aq such 
that for every 7 in Aq, there exist vertices 7', 7" adjacent to 7, with 
7 >- 7', 7", and such that 77', 77" have distinct slopes. 

For an ordering on Aq as in the previous lemma, we assign indices 
to the vertices in such a way that 7j >- 7^ when i > j. The first 
indices 1,2, . . . , /q are assigned to the vertices lying on the boundary. 
To those interior vertices which are joined to the boundary by two or 
more edges of distinct slope we assign the indice s f q + 1 n. Let us 



recall t he following notation and remarks from ISchenck and Stillman 



fll997a} ). For each interior vertex 7, and each / e JT{l), let /[7] denote 



the corresponding element in Hq{J). Then Hq{J') is generated by 
{/W I / ^ 1^(7) }• By definition of Ji^), we know that 

(13) /[7] = EC'^W 

for some polynomials G R. Notice that if r is an edge whose vertices 
are 7 and 7' then the 

(14) CV[7]=C[7'] 

in Hq{J)] in particular t^'^'^ f [7] = when r is an edge connecting 7 
to the boundary. 

Here is another lemma that we need to prove that iJo(i/fc) = for 
degree > 4r + 1: 

Lemma 5.2. Letii,i2 andi-^ be three equations of distinct lines through 
a point p and L the equation of a line that does not contain the point 
p. Then for any polynomial g of degree d > r + 1, there exist u,v & R 
of degree d and w & R of degree r — 1 such that 

g = u + e^^ V + L'^-^'+i r^' w, 

with w = in the case r = 0. 
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Proof. The case r = is direct from the hnear dependency of £i, £2 and 
£3. Suppose that r > 0. Since ii, £2 and L are hnearly independent, 
we can make the following change of coordinates: 

£1 = X] £3 = X + ay for some a 7^ ; 

£2 = y\ L = z. 

Let g = x'^y^z^ be a monomial such that i + j + k = d. li k < d — r 
then i + j >r and the polynomial £2^^g is in the ideal G {x^~^^ , y^~^^) . 
l{k>d — r + 1 then g is a multiple of z'^~^^^ = L'^~^^^. Thus we can 
write 

£1+^ g = £1+^ u' + £1+^ v' + L-^-'+i w', 

for some polynomials u',v' G i? of degree d and G -R of degree 2 r. As 
l^d-r+i^/ -g ideal £2'''"'^5 ^3'''"'^) and as L is a non-zero divisor 

modulo this ideal, we deduce that w' = £1+^ ^" + £3"^^ ^" + £3+^ ti?", with 
u",v'\w" G -R of degree r — 1. Collecting the coefficients of £'['^^,£2'^^, 
we obtain the desired decomposition: 

£-+1 g = £1+^ (u' + L'^-'-^^ u") + £^+1 [v' + L'^-^^+i t;") + L'^^^'+i £^+1 

We use this lemma to prove the following result: 

Theorem 5.3. The dimension of Cl{A) when k > Ar + 1 is given by 
the lower bound formula of Theorem \2.1[ 

Proof. From ([2]), and the formulas for the dimension of the modules 
0, Q and dZD, it suffices to show that HQ{J)k = for > 4r + 1. 
Equivalently, we need to show that f[^] =0 in Hq^J)^. for all 7 and 
/ G J'{'j)k- Ordering the vertices as in Lemma 15.1^ we consider the 
first interior vertex in the ordering, and denote it by 7. Let Ti,T2, . . . ,Tt 
be the edges (not necessarily with different slopes) containing 7 and 
coi, C02, cot be respectively the end points of Xj, which are distinct 
from 7. We number first the edges ri and T2 that connect 7 to the 
boundary, and the remaining edges counterclockwise, starting from T2. 
Let us denote by £±,£2, ■ ■ ■ ,£t (respectively) equations of these edges 
and by Li equations of the edges connecting cOi to Wj.i for i = 2, . . .t. 

We are going to prove by induction on j > 2 that for any homoge- 
neous polynomial / of degree > 3r, we have 

in HoiJ). 
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Figure 1 . Notation in Theorem [531 

Let us prove it first for j = 3. Let f & with A; > 3r. By construc- 
tion is on the boundary, thus we have L^'^^f [u^] = L^^^ f [(^2] = 0. 
By Lemma [5.21 

for some polynomials u,v,w E R. Then we have 

= Ll^+'r+'w[u,] = 0. 

Let us take now i > 2 and assume that the induction hypothesis is 
true for 3 < j < i. We consider first i^^^ f [uji] with / homogeneous of 
degree > 3r. By Lemma [5.21 applied to ii,i2,£3 and Lj, we have 

for some polynomials u,v,w & R. Then we have 

As L^'+' is in the ideal generated by i^^l, L^^l, by the induction 
hypothesis, we have 

This proves our claim for / [ui] = 0. 

Let us consider now L^^-^ f [coi] = L^^^ f [coi-i] with / homogeneous 
of degree > 3r. By Lemma [5.2[ 
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for some polynomials u,v,w G R. We deduce that 



fr+l JJ'+^ 







by the induction hypothesis and the previous computation. 

This concludes the induction proof and shows that for any / G 
jr(7)4r+i, we have 

1=1 i=3 

= 0. 

Therefore Hq{J) = and the dimension of C^(A) is given by the lower 
bound when k < 4r + 1. 

6. Examples and remarks 

In this section we compute the dimension formula for some com- 
plexes. We prove that the dimension formula is in fact given by the 
corresponding lower bound given in Theorem 12.11 In the first example 
we show directly that Ho{J')k = 0. In the second example we number 
the vertices in such a way that the upper bound we obtain agrees with 
the lower one. The subdiv isions we consider are tw o of the so-called 
Powell-Sabin subdivisions ( Powell and Sabin 



Example 6.1. Powell-Sabin 12-split. 

Let A be the simplicial complex supported on a triangle |A| subdi- 
vided into twelve smaller triangles, as in Figure [2l 



A / 






72 \ 




<7o / \ 




TV''* \ 











Figure 2. Powell-Sabin 12-split. 

Let us note that by using any numbering on Ag, the upper bound 
we get from Theorem 13.11 equals the lower bound, leading easily to the 
dimension formula. On the other hand, it is not possible to find a 
numbering on the interior vertices of A with the condition required to 
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apply Theorem 14.11 Hence it is not possible to find an upper bound 
for this spline space by using the formula of Schumaker. 

By the remarks in Section [5l it is easy to check that i^j'^^f [ji] = 
for every 7j and t^^^ f G J'i'Ji). Thus, for any k and any r < k, the 
homology module ifo(i7)fc is zero and the dimension of the spline space 
C^(A) is given by the formula in Theorem 12.11 

This example is an instance of a tri angulation where all edges ar e 
pseudoboundaries (in the terminology of lSchenck and Stillman fll997bh) 



and the computation also follows by Lemma 2.5 of lSchenck and Stillman 



U997M- 



Example 6.2. Powell-Sabin 6-split. 

Let A be the simplicial complex supported on a simply connected 
triangulated region |A| in M^. Assume the triangulation consists of 
/o vertices (interior and exterior) and /° i-dimensional interior faces 
(z = 0,1,2), where /2 is taken as the number of all the triangles cTj 
in the subdivision. We refine this triangulation by subdividing each 
triangle into 6 triangles to get A, as follows. 

(1) In each triangle G A2 choose an interior point i^i, in such a 
way that if two triangles (Tj, aj have an edge r in common, then 
the line joining i/j and uj intersects r at an interior point fiij. 

(2) Join each new point z/j to the vertices of the triangle Uj, and to 
the points fiij on the edges of cxj. 

(3) For the triangles having an edge (or more) on the boundary, 
choose a point in the interior of such edge and join it with i/j. 




Figure 3. Powell-Sabin 6-split of A. 
We want to apply Corollary 13.21 and give a formula for dimC^(A). 
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Let US consider the following numbering on the vertices in Aq. We 
take first a triangle cTj G A with (at least) one of its edges on the 
boundary, denote by 7^ the vertex in cxj that is in the interior of A. 
With the notation as above, we assign the index 1 to the vertex z/j, the 
indices 2 and 3 to the vertices /ijj on the edges of cxj, and finally the 
index 4 to the vertex 7j. After this, we consider the edges of cTj as if 
they were all part of the boundary and iterate the previous process. We 
choose another triangle in A having (at least) one edge on the "new" 
boundary. We index its vertices following the same order as we did for 
the ones in Uj. We continue this numbering until we have considered all 
the triangles in A and hence have index all the vertices in Aq. Figure 
m shows the triangles that could be considered first, and the different 
colors refer to the values of ti for the interior vertices in each case. 




Figure 4. Numbering on the vertices in Aq. 



Let us notice that for this numbering, the values for ti corresponding 
to these vertices are either > 3, as is the case of the vertices Ui and 7j, 
or ii = ti = 2, as it happens for the vertices fiij. Thus, we arrive to the 
following result. 



k~l\ f2k - 1\ , , 
2+2 l(/.-2). 



Proposition 6.3. For any k > 3: 

dimC,^(A)=(^ + 2)+/°(^2^)+2/Q° 
For k = 2: 

dimC2^(A)=3/o. 

Proof. It follows directly from the remarks we made above about tj 
for the counting, and Corollary 13.21 

For the numbering we defined on Aq, it is not possible to apply 
Theorem 14.11 to find an exact value for the dimension. In fact, when 
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the initial subdivision consists of more than one triangle, by using the 
numberings allowed in that theorem, we get upper bounds which are 
strictly bigger that the actual dimension of the space. 

Remark 1. In this paper we have confined our attention to triangular 
partitions but it is easy to check that our proofs of the lower and the 
upper bounds can be extended to rectilinear partitions. The construc- 
tion also allows to consider regions subdivided by curved boundaries, 
but then the ideals must be considered individually. The boundary of 
|A| itself can, of course, be curved. 

Remark 2. The resolution f or ideals generated by power of linear 
forms presented in the paper by Geramita and Schenck ( 19981 ). applied 



also when the power of the linear forms is different. It makes possible 
to apply the ideas we present in this paper to mixed splines, i.e., splines 
where the order of smoothness may differ on the various edges of the 
subdivision. 

Remark 3. The proof of the exact dimension formula when the degree 
of the spline space is > 4r + 1 cannot be directly extended to spaces 
of polynomials of degree k > 3r + 1. With some restrictions associated 
to the number of different slopes of the edges containing each vertex it 
is possible to prove the result for this degree. 

Remark 4. We hope that the two methods we used in the examples 
to find the exact dimension, namely, by showing that the lower and 
the upper bounds coincide for certain numbering on the vertices, or 
that Hq{J) = directly by considering the equations of the edges, 
illustrate the way to easily find the dimension of the spline space for 
many particular triangulations. 



Rema rk 5 . There is a version of th e upper bound given by ISchumaker 



( 1l984j ) in iLai and Schumakerl (120071 ) , where the constraint on the num 



bering of the vertices is relaxed. However, for this new version of the 
theorem, there are still examples for which the upper bound formula 
does not give the correct r esult, see for instance the example p. 243 of 



Lai and Schumaken ( 120071 ). 



Remark 6. For the three dimensional version of the problem, it is 
necessary to analyse ideals generated by powers of linear forms in three 
variables. The formula of the dimension, analogous to formula ([2]) in 
Section 121 have two homology modules, and in order to approximate the 
dimensions the dimension of these two modules must be bound. Some 
issues and positive results about this problem will be in an upcoming 
paper. 
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